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Using Schrodinger functional techniques, we determine the coefficient of the clover term nec- 
essary for non-perturbative O(a) improvement of hypercubic smeared Wilson fermions on a 
quenched plaquette action background. Unlike for unsmeared Wilson fermions, the resulting 
clover coefficients are close to the tree-level value even at coarse lattice spacings, indicating the 
absence of large cutoff effects. A number of exploratory tests are also performed with the im- 
proved action. 



X 



The XXV International Symposium on Lattice Field Theory 
July 30 - August 4 2007 
Regensburg, Germany 



* Speaker. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http://pos.sissa.it/ 



Non-perturbative improvement ofnHYP smeared Wilson fermions 



Roland Hoffmann 



1. Motivation 

We present results from a non-perturbative determination of the clover coefficient csw for 
Wilson fermions that couple to hypercubic smeared gauge links W. This is done on a quenched 
gauge background and while the actual results will have little practical relevance, they provide an 
interesting comparison to the unsmeared case. In this way we want to quantify how much smearing 
can help with Wilson-type fermions at coarse lattice spacings, especially concerning their chiral 
properties. 

We use the same gauge background (plaquette gauge action) and technique as the ALPHA col- 
laboration in their determination of c$w for thin links [[3p. The quantitative comparison of smeared 
and unsmeared case includes the value of the clover coefficient itself, remnant cutoff effects, the 
finite renormalization constant of the vector current and the range of accessible quark masses at a 
given lattice spacing. 

2. Setup and strategy 

Simulations of the quenched QCD Schrodinger functional (SF) are performed using the pla- 
quette gauge action and two quenched flavors of improved Wilson fermions. For unexplained no- 
tation regarding the hypercubic (nHYP) smearing or the Schrodinger functional we refer the reader 
to Refs. The SF setup imposes Dirichlet boundary conditions in the temporal directions and 

thus provides another IR cutoff in addition to the quark mass. Therefore, in principle, simulations 
at the critical point are possible. Moreover, the Schrodinger functional is a convenient formalism 
to formulate and implement non-perturbative improvement and renormalization conditions. 

The nHYP smearing under consideration here was recently used in dynamical Wilson clover 
simulations |§, ^] as well as in the construction of the kernel for both quenched and dynamical 
overlap studies 

When used with Schrodinger functional boundary conditions, the HYP construction [JIJ ||] 
has to be modified in the vicinity of the temporal boundaries. We adopt the prescription that the 
Dirichlet conditions are preserved under the blocking such that the original nHYP construction 
(labeling lattice sites by n) 

V HJI =Proj t/(3) [(l-a 1 )t/ n# + ^ £ Vn^Vn+^vl ] , (2.1) 

V n ,n;v = P r °jj/(3)[(1 ~ &2,)Un,il + ~T £ V n.py \N n+p ,\i;p vV n+fl,p;V p\ > ( 2 - 2 ) 

F„ #;vp =Proj t/( 3 ) [(l-a 3 )t/„ # + Y £ u n,n u n+n,nUl +{ln ]. (2.3) 

±r)^py,p. 

is supplemented by the prescription (labeling lattice sites by (xo,x)) 

V,(0,x) = Vi-(0,x) = y,(0,x) = Ui(0,x) (2.4) 
Vi(T,x) = Vi(T,x) = Vi(T,x) = Ui(T,x) (2.5) 

for the spatial links in the temporal boundaries. Here £/,-(0,x) and Ui(T,x) are the prescribed spatial 
boundary fields and we use (X\ = 0.75, (X2 = 0.6, 0:3 = 0.3 as blocking parameters. We note in 



2 



Non-perturbative improvement ofnHYP smeared Wilson fermions 



Roland Hoffmann 



passing that the classical background field induced by abelian homogeneous boundary conditions 
is invariant under (iterated) smearing as long as a projection/normalization is employed. 

The improvement condition itself consists of demanding that AM, the difference between two 
definitions of a current quark mass, is equal to its tree-level value AM ^ [0]. The cutoff effect AM 
is evaluated at the critical point (defined through yet another current quark mass) and the clover 
coefficient is tuned until this condition is satisfied. Note that the use of a color background field in 
the SF allows quark mass definitions that are 0(a)-improved without the knowledge of the axial 
current improvement coefficient ca [ff]. To summarize, at each value of the inverse bare gauge 
coupling j6 a number of csw values are simulated, where in turn for each of those the hopping 
parameter has to be tuned to find the critical point. An interpolation in csw then gives the point 
where AM = AM(°). 

The SF boundary improvement terms are dropped since they do not affect the observables we 
will be interested in Lastly, physical units for the results are obtained using the hadronic scale 
ro ffO| Ol anc ^ assuming tq = 0.5 fm. 



3. Numerical results 

All simulations are performed on lattices of size 8 3 x 16 and the background field is chosen 
as in [^J. First, we test the sensitivity of the improvement condition, i.e. how AM varies with 
the clover coefficient, with results from j3 = 6.45 shown as filled symbols in Fig. |]. The point 
csw = 1-043(11), where AM is equal to its tree-level value 0.00028, is also indicated. We obtain a 
clear signal for the clover coefficient with a value that is rather close to one. 

For comparison we also show results from [|3|] at a similar lattice spacing (j8 = 6.4, open 
symbols). Clearly, the sensitivity is very similar, while the value of the clover coefficient, that 
is necessary to achieve non-perturbative improvement, differs greatly. In other words, the cutoff 
effects that are being canceled with a non-perturbative csw are dramatically reduced if the fermions 
are coupled to the nHYP smeared gauge field. 




c sw 9q 



Figure 1: Sensitivity of the improvement condi- Figure 2: Non-perturbative result for csw with 
tion: AM vs. csw for nHYP links at /3 = 6.45 nHYP links and comparison to the thin link re- 
plied symbols) and thin links at j3 = 6.4 (open suits, 
symbols, [^J). 
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The determination was done for 5.8 < /? < 12 and the resulting non-perturbative values for 
csw are plotted in Fig. ||[ Note that, while the error on our determination of the clover coefficient 
grows when one goes to coarser lattices due to the increased statistical fluctuations, the value itself 
rises only slowly. At j3 = 6 we obtain csw = 1.107(15) compared to c sw ~ 1.8 with thin links |^]. 
Our data is well described by the interpolating formula 



csw(P) 



p -5.611 
j8 -5.647 



(3.1) 



By design the result from the non-perturbative improvement condition will approach the per- 
turbative expansion for small enough value of gp. From the plot it is then evident that the 1-loop 
contribution to csw is very small in the case of nHYP smearing. More precisely, from the asymp- 
totic behavior of Eq. ( |3~T| ) we would estimate ~ 0.006 as opposed to = 0.266 without 
smearing. 

After fixing csw non-perturbatively Ref. @ proceeds with the improvement of the axial cur- 
rent, i.e. a determination of the coefficient ca- Already with thin links the remaining cutoff effects 
are rather small and so is the value of ca- Various methods to determine the axial current improve- 



ment coefficient are available [ |12| ] and we have tested those that utilize variations in the quark mass 
with respect to the insertion time and the periodicity angle d of the spatial boundary conditions. 
No sensible criterion could be found since the observed cutoff effects were exceedingly small, 



but a more detailed study and/or the use of wave functions [ ]13| ] might lead to better results. We 
therefore proceed to test the quenched nHYP clover action with c$w given by Eq. (3. 1) and ca = 0. 
The following simulations are performed with vanishing background field. 



4. Testing the improved action 

4.1 Remnant cutoff effects 

As a first test we consider the time dependence of the PC AC quark mass on a 16 3 x 32 lattice 
at j6 = 6.2, corresponding to a lattice spacing of a = 0.068 fm. In Fig. || the result is compared to 
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Figure 3: The PCAC mass at j3 = 6.2 with and 
without nHYP smearing after non-perturbative 
improvement. 




Figure 4: The PCAC mass at j3 =6 with nHYP 
smearing for two values of the spatial fermionic 
boundary phase. 
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a thin link simulation from Ref. [[J at a similar mass. In both cases the quark mass is contained 
in a narrow band of ~ 4MeV width in the interior of the lattice. However, toward the temporal 
boundaries the smeared action clearly shows smaller cutoff effects since it remains much natter 
than in the thin link case. 

Since for a smeared action this is an already rather fine lattice, we also perform a test at 
j8 = 6, a ~ 0.1 fm, and 8 3 x 16. Fig. || shows the PC AC quark mass with non-perturbative csw at 
fc = 0.1257 and 6=0 and 1. By themselves, both data sets are contained in a 3MeV band in the 
interior of the lattice and even together, the spread is only ~ 5 MeV. For reference, note that the 
mass difference between 6=0 and 1 at xq = T /2 in lattice units is only 0.00079 or roughly twice 
its tree-level value. 

4.2 Vector current renormalization 

Let us now consider the renormalization of the (local) isovector vector current. In a mass- 
independent renormalization scheme the latter is given by 



(V R )^=Z v (l+byam q )V£ 



(4.1) 



where V£ is the bare current, m q the bare subtracted quark mass and the term proportional to 
by is required to preserve 0(a)-improvement under renormalization [f|]. Both Zy and by have a 
perturbative expansion of the form 1 +0(gg). For the thin link case the finite and scale-independent 



renormalization constant Zy was calculated in [ 14] and found to deviate significantly from its tree 



level value at relevant lattice spacings. The same also holds for the dynamical case [15]. 

Zy itself is obtained from a ratio of SF correlation functions at zero quark mass, while by can 
be calculated from the slope of this ratio as a function of the subtracted quark mass m q at m q ~ 0. 
Results are shown as the filled symbols in Fig. |5[ Note that these were not obtained on a line of 
constant physics but at fixed L/a and thus include a (presumably small) 0(a 2 /L 2 ) uncertainty. In 
both cases the results using nHYP smearing are significantly closer to unity than with thin links. 
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Figure 5: Non-perturbative results for Zy and b\ for nHYP smeared Wilson fermions (filled circles, error 
bars are smaller than symbol size). For comparison the interpolating formulae for thin links from Ref. [ |l4| | 
are shown as dashed lines. 
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Even at a lattice spacing of 0. 1 fm (gg 
current is almost conserved. 



1) we have Zy ^ 0.97, indicating that also the local vector 



4.3 Exceptional configurations 



Finally, we come to the issue of exceptional configurations. A tree-level analysis [16] has 
shown that the Schrodinger functional boundary conditions induce an IR cutoff proportional to 
1/r 2 in the spectrum of the squared Hermitian Dirac operator. But due to the absence of chiral 
symmetry, the Dirac operator of the interacting theory can (and does) develop eigenvalues much 
smaller than this bound. These lead to exceptional configurations where estimators for fermionic 
correlation functions have values that are orders of magnitude above the normal level of fluctua- 
tions. 

In practice it was found [^] that in the quenched SF with plaquette gauge action /3 > 6.4 is 
necessary to avoid exceptional configurations at zero quark mass. Their occurrence is known to be 
linked to extremely localized fluctuations of the gauge fields, so called dislocations. The Wilson 
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operator, especially with a clover term, is sensitive to those [17] and since the hypercubic smearing 
was designed to optimally suppress dislocations, it should be well suited to address this issue with 
Wilson clover fermions. 

In Fig. |6| we show the current quark 
mass as a function of (2k:) -1 from 8 3 x 16 
simulations at j8 = 6, 6 = and vanishing 
background field. No sign of exceptional 
configurations was found in the statisti- 
cal analysis. Repeating this test on even 
coarser lattices (/3 = 5.8, a = 0.136fm) 
shows that there a quark mass of about 
25MeV is required to avoid exceptional 
configurations. Strengthening the SF cut- 
off with a smaller time extension T ja = 
8 again allows simulations at the critical 
point even at this very coarse lattice. To 
summarize, we estimate that in the SF the 
lattice spacing that is accessible along the 
critical line is roughly doubled (j8 = 6 vs. 
6.4) when nHYP smearing is employed. 
We also note that at the same time the additive mass renormalization is reduced by almost an order 
of magnitude. At j8 = 6 we obtain fQ. = 0.12617, corresponding to am c = -0.037, whereas Ref. |g] 
quotes K c = 0.135196, i.e. am c = —0.302. 
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Figure 6: The current quark mass as a function of (2 k - ) -1 
with nHYP smearing at j3 = 6. The data points are labeled 
with the corresponding hopping parameter. 



5. Conclusion 



In this short study we have successfully used nHYP smeared clover fermions in the frame- 
work of the Schrodinger functional. We implemented the non-perturbative action improvement 
pioneered by the ALPHA collaboration and found that no large values of the clover coefficient are 
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required to cancel the cutoff effects under consideration. Instead, csw remains close to its tree-level 
value even on coarse lattices. 

As an example of a finite renormalization constant we performed the non-perturbative renor- 
malization of the local vector current. The renormalization constant Zy differs from unity by only 
3% at a lattice spacing a = 0.1 fm. 

We hope to have convinced the reader that also with modern (smeared) fermion formulations 
the Schrodinger functional remains a useful tool for non-perturbative improvement and renormal- 
ization. Once a dynamical implementation is available, it will be very interesting to check to what 
extent smearing helps with respect to the large remnant 0{a 2 ) effects and ambiguities found in two 



flavor clover simulations [15,18 
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